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Introduction 



Astrochemistry and molecular astrophysics are interdisciplinary fields involving chemistry, physics, and astronomy. Their main 
aim is to understand the chemical evolution of the universe. Astrochemistry and molecular astrophysics are often used as synonym 
and a distinct line between them is difficult to place. Both research areas study the formation and destruction of molecules in space, 
their interaction with radiation and their feedback on physics of the environments. As a sort of simplification, it can be affirmed 
that astrochemistry focuses on the study of the chemical processes at work in space and on the understanding of the molecular 
evolution, trying to address the question of molecular complexity. On the other hand, molecular astrophysics is more interested in 
the information that can be retrieved from the study of molecules in the universe. In fact, molecules being ubiquitous (they are 
found everywhere in space: in interstellar, circumstellar and pregalactic gas, in protostellar disks, in the atmospheres of stars and 
planets, . . .) are unique probes to retrieve information on molecular excitation, radiative transfer, and kinematics. In a few words, 
this simplification affirms that astrochemistry studies the evolution of molecules abundance, while molecular astrophysics studies 
the molecules as physical diagnostics. In other words, we can say that astrochemistry is the study of the abundance and reactions of 
chemical elements and molecules in the universe, and their interaction with radiation, while molecular astrophysics concerns the 
study of emission from molecules in space. The line between astrochemistry and molecular astrophysics being so fine and vague, 
will not be addressed further, while the two fields will be considered as parts of the same research areas. 

The approach at the basis of astrochemistry and molecular astrophysics is depicted in Figure 1. Both fields span astronomical 
observations, modeling, and theoretical and/or experimental laboratory-based investigations. From this figure it is apparent how 
wide the corresponding research field is; consequently, it cannot be exhaustively addressed in this contribution. Interested readers 
are therefore referred to the 'Further Reading' section, whereas we limit our discussion to the role played by quantum chemistry in 
the field of astrochemistry as well as to the contribution of experimental rotational spectroscopy to molecular astrophysics. 
Quantum-chemical calculations are mainly employed in the 'laboratory' studies: spectroscopic characterizations that may support 
astronomical observations as well as investigations of chemical reactions (also including processes on grains) to determine 
thermochemistry, kinetics, and reaction mechanisms. In the latter context, actual studies are those that aim at understanding the 
level of complexity that can be reached by molecular species in space. Experimental laboratory studies range from the field of 
molecular spectroscopy (which involve observation of transitions between different energetic levels, either rotational or vibrational 
or electronic) to simulations of the reaction conditions in extraterrestrial environments (which involve, e.g., the study of 
ion-neutral reactions, of reactions taking place in either aerosols or cold plasma, etc.). 

About 180 molecules have been detected in space, with at least 175 molecules observed unambiguously in the interstellar 
medium (ISM) or circumstellar shells and more than 50 molecules detected in extragalactic sources (see, e.g., the Cologne 
database 1,2 ). Our knowledge on the universe chemical inventory has been obtained and continuously updated by means of 
astronomical observations, which have revealed a plethora of molecular species, 3 thus demonstrating the rich chemistry of the 
universe. Molecules can exist in a wide range of astrophysical environments, from the extremely cold regions between stars to the 
atmospheres of stars themselves. Gas-phase species have been mostly discovered via their rotational signatures, with the corre- 
sponding frequencies going from the millimeter-wave region up to the far-infrared, even though spectroscopies in the visible and 
infrared domains have also led to the detection of molecular species. The observation of gas-phase species performed using 
rotational spectroscopy can be undertaken from the ground, 4 while vibrational spectroscopy observations, for which ground-based 
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Figure 1 The approach at the basis of astrochemistry and molecular astrophysics. 

observations are hampered by the terrestrial atmosphere, are mainly carried out by spectrometers on board satellites or airbornes. 
In particular, infrared spectroscopy is suitable for retrieving the chemical composition of planetary atmospheres. As concerns the 
dimension, the molecules detected ranges from diatomics to species containing more than ten atoms, up to polycyclic aromatic 
hydrocarbons that were identified to be abundant, for instance, in the atmosphere of Titan, Saturn's largest moon. The detection 
also includes isotopic species. Up to now, most of the molecules that have been detected in the ISM or elsewhere in space are 
neutral species. Some cations have also been found, along with only a few anions whose discovery is quite recent. 

This contribution is organized as follows. In the next section, a brief overview on laboratory studies with emphasis on 
computational chemistry and molecular spectroscopy is provided. Thereafter, selected exemplificative cases are presented and 
discussed. The contribution is completed with a 'Further Reading' section. 



Laboratory Studies 

Computational Chemistry 

As mentioned above, in the field of astrochemistry and molecular astrophysics quantum-chemical calculations can be used to 
provide spectroscopic characterizations that can guide the observation in the laboratory or even directly in space. To support 
experimental investigation, theoretical predictions should be accurate and thus imply high-level quantum-chemical calculations to 
be used. The focus is here on rotational and vibrational spectroscopies. The theoretical prediction of the corresponding spectra 
involves two tasks, namely (a) the determination of the required spectroscopic parameters, and (b) the simulation of the spectra 
based on the given set of spectroscopic parameters. The latter typically involves the diagonalization of an appropriate Hamiltonian, 
while the accurate determination of the spectroscopic parameters requires a quantum-mechanical treatment of the molecular 
system under consideration. To fulfill the accuracy requirements, the key point is the employment of the coupled-cluster (CC) level 
of theory. In particular, the CC singles and doubles (CCSD) approximation augmented by a perturbative treatment of triple 
excitations (CCSD(T)) 5 provides a good compromise between accuracy and computational cost. On the other hand, for large 
systems methods rooted in the density functional theory are demonstrated to provide results that meet the spectroscopic accuracy 
requirements. 6 For exhaustive accounts on the quantum-chemical calculations of rotational parameters, the readers are referred to 
Refs. [7-9], while for vibrational spectroscopy Refs. [10-15] (and references therein) provide an overview of methodological and 
computational approaches (for both perturbative and variational models). Here, we only briefly mention that composite schemes 
are largely employed in order to account for various contributions (basis-set and electron-correlation effects) separately at the 
highest possible level and to combine them for obtaining the best possible theoretical estimates (see, e.g., Ref. [16] and references 
therein). Similar computational requirements are needed for accurate thermochemistry investigations: CC techniques in conjunc- 
tion with composite schemes permit to obtain a sub-kj mol -1 accuracy (see, e.g., Refs. [17-20]). 

For an accurate modeling of the chemistry that takes place in chemistry space, the understanding of chemical kinetics in 
astrophysical environments is a challenge and plays an important role. In fact, it requires to predict reaction rate coefficients for 
temperatures (250-2000 K) and conditions well outside those commonly encountered in terrestrial laboratories. Quantum-chemical 
calculations can also be employed in the field of chemical kinetics, not only in the modeling of the potential energy surface (PES) (thus 
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pointing out reaction barriers to be overcome) but also in the determination of reaction rate constants. For the latter, stochastic or 
(semi-) classical approaches can be employed. On general grounds, master equations should be formulated and then solved either 
analytically or with numerical techniques. Refs. [21,22] provide explicative examples involving systems of interest to the present 
context (the HO + CO reaction and the formation of glycine, respectively) . 

Molecular Spectroscopy 

As mentioned above, molecular species detected in space have been mostly identified via their rotational signatures, which in turn 
were previously measured in laboratory or obtained from the rotational parameters determined from laboratory studies. Due to its 
relevance in the field of molecular astrophysics and the limited account here reported, only rotational spectroscopy will be 
considered in this contribution from the experimental point of view. 

Rotational transitions typically lie in the range of 3 GHz-3 THz, 25 with rotational constants being the parameters that mainly 
govern a rotational spectrum. Therefore, detections in the millimeter and submillimeter wavelength range are usually limited to the 
study of small- to medium-sized molecules, while measurements in the centimeter-wave region allow one to investigate larger 
molecules, even of biological relevance, and molecular complexes. Rotational spectroscopy can be performed either in the time or 
frequency domain. In the first case, the typical resolution and frequency accuracy are about 5-20 kHz and ~1 — 5 kHz, respectively, 
while in the latter case the resolution is mainly limited by Doppler effect, the resolving power being 50 kHz-2 MHz, according to 
the frequency range considered. The Lamb-dip technique 24 is a powerful tool for obtaining sub-Doppler resolution in the high 
frequency regime, even up to the terahertz frequency range. 25 

Examples and Discussion 

As mentioned above, too vast are the fields of computational astrochemistry and molecular astrophysics to provide an exhaustive 
overview in this contribution and several are the possible applications of quantum-chemistry and/or molecular spectroscopy. 
In the following, an attempt of providing a brief but significative collection of examples is made. First, computational studies will 
be addressed, while in the second part of this section examples of applications of rotational spectroscopy to the field of Astrophysics 
will be discussed. 

Computational Studies 

Computational chemistry allows us to infer how complex molecules can be formed in the ISM, to understand which molecular 
species exist in space (planets' atmospheres, ISM, etc.), and to support or to provide the spectroscopic characterization required for 
astronomical observations. 

Glycine: a potential interstellar biomolecule 

The existence of biological building-block molecules in the ISM has attracted a great interest, especially because of their direct 
connection with the chemical evolution of the universe. Glycine (NH 2 CH 2 COOH) is the simplest amino acid and therefore it has 
been intensively searched for in the ISM. While glycine has been discovered in some meteorites fallen in Earth (e.g., Murchison 
meteorite 26 ) and within the returned samples of comet 81P/Wild2 from the Stardust mission, 2 ' its detection in the ISM has been 
debated. Despite several experimental works, direct measurements to elucidate reaction mechanisms are difficult to perform, while 
quantum-chemistry can be used to shed light on the mechanistic steps of the amino acids formation in the ISM by exploring the 
corresponding PESs. 

Recent studies (see Ref. [16] and references therein) demonstrated that state-of-the-art quantum-mechanical approaches permit 
to accurately and fully characterize the conformational behavior of glycine, with conformational enthalpies determined well within 
the chemical accuracy. As is evident from Table 1 of Ref. [16], the gas-phase thermodynamic properties were obtained with a 
1 kj moP 1 accuracy. It is worthwhile noting that the same accuracy also applies to the transition state structures connecting the 
local minima on the PES. Such result is an important step toward the identification of the formation mechanism of glycine in 
the ISM. 

Rimola et al. 22 recently proposed a possible path for the formation of glycine-derivative species which may explain the discover 
of glycine in some meteorites and comets. The water-dominated ice mantle covering a silicate-based core of an interstellar grain was 
simulated by means of a cluster of eight water molecules. The key point of the reaction mechanism is the formation of OH surface 
radical species due to ionizing particles, such as cosmic rays. If ice contains CO and methylenimine (NHCH 2 ), the reaction occurs 
via the formation of the COOH radical due to the coupling of the OH radical with an incoming CO molecule. Subsequently, by 
reacting with NHCH 2 the COOH radical yields a radical glycine-derivative species, which finally leads to glycine by the addition of 
one H atom. While interested readers are referred to Ref. [16] for a complete account, here we point out that quantum-chemical 
calculations, in addition to explain the formation of glycine on interstellar grains, can also provide information on the reaction 
rates. In fact, as pointed out above, the Gibbs free energy of the reaction transition state required in Eyring's equation can be 
accurately determined, thus allowing one to estimate the reaction rate at a given temperature with good accuracy. This is indeed a 
key information in order to indicate whether reactions are feasible or not on astronomical time scales. 
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Organic molecules in Titan's atmosphere: cyclopropenyl cation and oxirane 

There is a great interest in Titan's atmosphere due to the fact that recent measurements revealed that Titan's atmosphere is 
characterized by a rich and complex organic chemistry (see Ref. [28] and references therein). Titan is thought to represent a 
model of primitive earth because its atmosphere contains significant quantities of carbon (^2% CH 4 ) and nitrogen (98% N 2 ), and 
traces of oxygen (~50 ppm CO). 29 Photochemical production of complex molecules containing C, N, O, and H is therefore 
possible, and consequently, the investigation of Titan's atmosphere might provide an unique opportunity to shed light on the 
organic evolution in the atmosphere of early Earth. 

The ion and neutral mass spectrometer (INMS) onboard the Cassini spacecraft revealed the presence of carbocations, such as 
CH 3 + , CH 5 + , C 2 H 3 + , C 2 H 5 + , C 3 H 5 + , and C 3 H 7 + , in Titan's atmosphere. In Ref. [28], based on INMS Cassini observations, on 
the analysis of laboratory crossed molecular-beam measurements, and on ab initio studies of the reactions of CH 3 + with ethylene 
and acetylene, it was inferred that a significant fraction of the C 3 H 3 + detected by INMS in Titan's upper atmosphere is the 
cyclopropenyl cation, which is the simplest Huckel's aromatic molecule. Furthermore, an important class of reactions in Titan's 
upper atmosphere, such as those of CH 3 + with methylacetylene and dimethylacetylene, were identified. The latter reaction explains 
the observed composition of C 4 H 5 + and C 5 H 7 + in terms of methyl-substituted and dimethyl-substituted cyclopropenyl cations, 
respectively. 

Huckel's cyclic molecules are stabilized by aromaticity. Formation of aromatic compounds is an important issue because the 
chemistry of Titan's upper atmosphere may shed light on the first steps toward a mechanistic theory of the terrestrial prebiotic 
chemistry. In fact, it has been largely accepted that Titan could represent a model of primitive earth. 29 Investigation of Titan's 
atmosphere requires quantum-chemical calculations of the energetics of the reactions and molecular species involved as well as 
astronomical (rotational as well as infrared) observations. The latter in turn are based on either experimental or theoretical 
laboratory studies. 

Recent photochemical models suggest that oxygen should be incorporated into the chemistry of Titan's atmosphere. 30 
Experimental simulations of Titan's atmosphere which includes CO at the 100 ppm level on N 2 -CH 4 mixtures led to the detection 
of an important prebiotic species, oxirane (C 2 H 4 0, also named ethylene oxide), in the products mixture. Since INMS is not able 
to mass resolve the CO molecule from the major species N 2 and, on a similar basis, heavier molecules or ions, the identification of 
oxirane in Titan's atmosphere is demanded to astronomical observations either in the infrared region or in the millimeter/ 
submillimeter-wave regime. The required spectroscopic characterization can be carried out either experimentally or theoretically. 
Focusing on rotational spectroscopy, Figure 2 demonstrates the accuracy that can be reached nowadays by high-level quantum- 
chemical calculations: for selected rotational transitions, experiment 12 and theory 33 are compared. It is noted that the calculated 
transitions are on average underestimated by about 0.1%. Therefore, theoretical predictions are sufficiently accurate for either 
laboratory or astronomical assignments. 

Moving to the field of vibrational spectroscopy, once again oxirane offers an interesting example. Bernstein and Lynch 
demonstrated a close correlation between the IR spectra of oxirane and cyclopropenylidene (c-C 3 H 2 ) with the major unidentified 
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Figure 3 The comparison in the 2.4-4.8 \im wavelength range between the simulated emission spectra of oxirane (red, lower trace) and the observed 
spectra (black, upper trace; from the uniform database of 2.4-45.4 urn ISO-SW spectrum of planetary nebula NGC 7027 35 ) is shown. Identified line 
fluxes reported by Salas et al. 36 are marked by gray lines. 



infrared (UIR) bands. 34 However, they also pointed out considerable uncertainties in some of the assignments and therefore called 
for a quantitative simulations of UIR spectra that properly account for anharmonicity effects on both line positions and intensities. 
Figure 3 shows the simulation of the IR spectrum of oxirane in the 2.4-4.8 um based on high-level quantum-chemical compu- 
tations, also accounting for anharmonic effects, 31 compared with the corresponding observation (ISO-SW spectrum of the 
planetary nebula NGC 70 2 7 35 ). The portion depicted also shows all identified line fluxes. 16 We note that the intense features of 
oxirane permit to explain the observed UIR features in the 3.2-3.4 (tm range. 



Astronomical detection of ions in the ISM 

The presence of negative ions in the ISM was predicted many years ago on general grounds or on the basis of ion-molecule chemical 
models (see, e.g., Ref. [37]). It was pointed out that a high electron affinity and a large number of vibrational states increase greatly 
the sticking coefficient of electrons; for instance, large, negatively charged carbon chains of the form C„H~ were predicted to be 
quite abundant. However, for a long time negative ions escaped detection because of the lack of accurate transition frequencies. 
Only in 2006, the first negative molecular ion C 6 H~ was detected in the laboratory and identified in the molecular envelope of IRC 
10216 and in the dense molecular cloud TMC-1. The spectroscopic constants derived from laboratory measurements were found 
identical to those derived from the astronomical data, thus leading to the unambiguous identification of C S H~ in the ISM. 

The missing observation of anions until 2006 was due to the lack of accurate transition frequencies. Quantum-chemical 
calculations played a crucial role in the identification of C 6 H~ in the ISM; in fact, the laboratory observation was guided by 
CCSD(T) computations. In 2008, a series of lines detected in the envelope of the C-rich star IRC+10216 was assigned to C 5 N~. 39 In 
this case, the assignment based on chemical and spectroscopic argumentations was confirmed by a perfect match with highly 
accurate computed rotational parameters. 



Rotational Spectroscopy 

As already mentioned, laboratory measurements of rotational transitions permitted the subsequent astronomical detection of most 
of the molecules discovered in space. In the following, a few exemplificative cases have been selected in order to present further or 
innovative applications of rotational spectroscopy in the field of molecular astrophysics. 



THz measurements 

Herschel Space Observatory, Stratospheric Observatory for Infrared Astronomy (SOFIA), and Atacama Large Millimeter Array 
(ALMA) opened up the submillimeter-wave region of the electromagnetic spectrum by making astronomical investigations of 
unprecedented high sensitivity and angular resolution possible. If on one hand the great spectral resolution that can be exploited by 
ALMA is expected to play a key role in the detection of complex species as well as rare isotopologues, on the other hand already the 
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first spectra recorded by Herschel turned out to be full of signatures, mostly unassigned. Unfortunately, a few molecular species 
show extremely numerous and relatively strong transitions throughout the submillimeter-wave region. Consequently, these species 
are likely to be a serious challenge as their emission occupies a substantial part of the spectral ranges available, thus impeding the 
study of other important species. The discovery of new species of astrophysical and astrobiological interest is therefore severely 
hampered unless the spectral lines from these 'weeds' can be removed from the data. For these reasons, accurate laboratory 
measurements in the submillimeter-wave and THz frequency ranges are required, and a great effort has been made in this direction. 

Transition frequencies in the region of the electromagnetic spectrum mentioned above require to be known with an accuracy 
preferably better than 100 kHz, thus leading to the need of extending laboratory data. Furthermore, extrapolations from 
low-frequency laboratory measurements might provide inaccurate and/or not reliable higher frequencies. An exemplificative 
case is provided by the experimental determination of the sextic centrifugal-distortion constant of HCO + (formyl ion), the most 
abundant molecular ion in the ISM, observed and observable in many objects. The results summarized in Table 2 of Ref. [40] point 
out the importance of accurate measurements in the submillimeter-wave region up to the THz in order to correctly determine the 
sextic centrifugal-distortion constant. In a previous work the attempt of evaluating it even provided a value with the wrong sign. The 
importance of a correct determination of these small spectroscopic terms is due to the fact that they play a key role in the prediction 
of rotational transitions at higher frequency, which in turn can be used for the identification in space of the molecular species under 
consideration. In passing, we note that the determination of unreliable spectroscopic parameters can be avoided if accurate 
computational data are used to support the fitting procedure. 

Accurate rest frequencies 

Although the available accuracy for most interstellar lines (~100 kHz) is sufficient for confident identification, detailed studies of 
molecular excitation, radiative transfer, systematic velocity gradients, and ambipolar diffusion in star-forming regions need 
significantly higher accuracy. Furthermore, the accuracy of astronomical observations is strongly influenced by the uncertainties 
of the various Doppler corrections due to the Earth motion. The best way to avoid such errors is to measure simultaneously the 
studied lines and the line of some other molecule with precisely known rest frequency (frequency standard). In Ref. [41], it was 
shown that rotational spectroscopy by means of Lamb-dip technique can fulfill the accuracy requirements for establishing 
frequency standard, because it permits to retrieve transition frequencies with an accuracy of about 1 kHz. In detail, in Ref. [41] 
the Lamb-dip technique allowed us to perform the most precise laboratory measurements of the rotational transitions of C ls O, 
thus leading to very accurate spectroscopic constants that permit to consider C ls O as a frequency standard for radio astronomical 
spectroscopy. While the main CO isotopologue cannot be used as a frequency standard because its rotational lines have extremely 
high opacity, C ls O was selected because its lines are usually optically thin. 

The discussions reported above and in the previous section point out the need of accurate measurements of rotational 
transitions, with particular emphasis on the THz region. While sub-Doppler resolution and the accuracy that this implies are 
well consolidated in the millimeter/submillimeter-wave region, completely different is the situation in the THz frequency range. 
Despite the explosive growth undergone by the corresponding technology in the last decade, the achievement of sub-Doppler 
resolution in this frequency region is only at its beginning (see Ref. [25] and references therein). In Refs. [25,42], it has been 
demonstrated that an accuracy of 1 kHz (i.e., an accuracy better than one part in 10 9 ) can be obtained in the THz region and how 
this accuracy permits to improve the spectroscopic parameters determination. 

Colli si onal broadening of rotational lines 

Collisional and radiative excitations are two important mechanisms which determine the existence and shape of spectral lines in 
the ISM. Therefore, collisional rates and transition moments are needed in order to interpret interstellar spectra in terms of local 
physical conditions. Experimental estimates for collisional rates rely on measurements of collision-induced spectral pressure 
broadening. 43 A significant example is provided by the investigation on N 2 H + perturbed by He, reported in Ref. [44], N 2 H + , one 
of the first molecular ions detected in interstellar space, plays an important role in understanding the physical conditions in 
protostellar cores, and has thus been the subject of extensive studies. From a computational point of view, the accurate knowledge 
of the PES for the N 2 H + -He interaction allows calculations of state-to-state collisional transitions rates, which are then used for 
probing the observations from interstellar molecular clouds. While a direct check of the accuracy of these calculations is not 
possible, the scattering matrix used can be validated by means of pressure-broadening and pressure-shift measurements. 44 These 
are measurements in which the collisional width and frequency are determined at different pressure values of the perturber (He in 
the example considered) and the corresponding pressure-broadening and pressure-shift coefficients derived. 
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